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We present two new approaches for the design of a Zeeman-Slower, which rely on optimal com- 
pliance with the adiabatic following condition and are applicable to a wide variety of systems. The 
first approach is an analytical one, based on the assumption that the noise in the system is position 
independent. When compared with the traditional approach, which requires constant deceleration, 
for a typical system, we show an improvement of ~ 10% in the maximal capture velocity, allowing 
for a larger slower acceptance, or a reduction of sim 25% in slower length, allowing for a simpler 
design and a better collimated beam. The second approach relies on an optimization of a system 
in which the magnetic field and the noise profile are well known. As an example, we use our 12-coil 
modular design and show an improvement of sim 9% in maximal capture velocity or, alternatively, 
a reduction of ~ 33% in slower length as compared with the traditional approach. 

PACS numbers: 07.55.-w,03.75.Be,39.10.+j,39.90.+d 



Systems utilizing atomic beams frequently require those beams to be decelerated, for example, in order to reduce 
the beam velocity to below the capture velocity of a magneto optical traps. Several schemes for slowing down 
neutral atomic beams are known, these include mechanical slowing collisions with cold background gas 0, pulsed 
laser fields [3|, pulsed electric fields [4], pulsed magnetic fields [5|, and Zeeman deceleration [6|, the subject of this 
communication. The Zeeman slower operates by matching a spatially varying magnetic field to compensate for the 
change in the doppler shift of the decelerated beam, permitting the use of a single frequency laser. 

The profile of this magnetic field is usually selected such that the atoms undergo constant deceleration. A design 
for an optimal coil shape which will produce this field has been presented @, but to the best of our knowledge, there 
is no published information examining the optimum field shape. Even though the constant-deceleration approach is 
mathematically simple, we show that it is far from ideal, and that other approaches are capable of optimizing various 
parameters of the experimental system. We present here two new approaches, both of which rely on satisfying the 
'Adiabatic Following' [8] condition in the best possible manner. 

The first approach is an analytical one, similar to the commonly used constant deceleration approach and is useful 
for the design of systems where the noise is independent of position along the slower. 

The second approach allows to optimize an already designed system where the magnetic field profile and noise are 
well known. As an example of such a system we consider a design used in our lab which aims to slow a metastable 
Neon beam for the purpose of trapping it in a magneto-optical-trap. 



The dissipative force exerted by classical light with wave- number k, on a two level system with line- width T, in a 
steady state is 



where sq — h/Is is the saturation parameter, Ii is the laser intensity, I s = hcTk 3 /12ir is the saturation intensity, and 
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is the doppler and Zeeman shifted laser detuning, where So is the difference between the laser and the atomic transition 
frequency and fi' is the magnetic moment of the transition. 

By requiring the detuned laser to be close to atomic resonance, S ~ in <j2j) we arrive at an approximate relation 
between velocity and magnetic field, 



v(z)= 6 f + ^B(z), (3) 

In order to slow the atomic beam from a initial velocity Vi to some final velocity v/, taking an on resonance laser, 
5q = 0, we immediately obtain the extrema of the required magnetic field, 

Brnax = B (z = 0) = ^rVi] B min = B (z = Z max ) = ^jV f . (4) 

M I* 

Introducing an overall detuning of the laser Sq ^ is compensated by an overall shift of the field, this is used for 
example when designed a 'midfield zero' (also called a 'spin flip') slower which has some advantageous features @. 
Nevertheless, from hereon we assume So = 0, the general result is easily derived by the substitution B — > B — %-Sq. 

Eq. ^ related the magnetic field extrema to two important parameters, the maximal velocity, Vi, and the length 
of the slower, z max . 

It is apparent that Vi should be large in order to slow a large fraction of the initial velocity distribution of the atomic 
beam, the need for a short slower is a somewhat more subtle. When interacting with the laser the atoms experience 
a random walk in velocity space resulting from the spontaneous emission of photons from the laser. This results in a 
non zero RMS transverse velocity or 'transverse heating' that contributes to decollimation of the beam. We are thus 
driven to make the slower as short as possible. 



III. ADIABATIC FOLLOWING 



Eq. relates the deceleration to the field gradient, however, Eq. (JTJ) enforces a strict limit on the maximum 
deceleration, 

_ hkT Sq 

Q"max — ~~Z ~ j • [pj 

By imposing the condition, a — — < a max and substituting velocity with the field using ^ we obtain a 
constraint on the field, 



9B Ojynax /«\ 

~ hk ~dz ~ Wb' 

hk 

Equation ^ is the 'adiabatic following' condition, see Q for a more thorough discussion. 

The magnetic field profile used should span from B max to B m i n while maintaining the condition ([5]). It is tempting 
to equate the two sides of ©, and solve to get, 



hk 



B(z) = — yj vf - 2a rnax -z. (7) 
This solution is also derived by demanding a constant and maximal deceleration. The slower length is found by the 

2 _ 2 

minimal field from Q to be z max = ^- — -. 

While this allows a short slower, any fluctuation from this field will result in violation of ^ and loss of the atoms 
from the slowing process. In order to take into account such fluctuations we introduce a positive 'noise parameter', 
a, which parametrizes the stability of the system under such fluctuations, 

a[z)= o^ + £dB (g) 
£-B hk oz 

hk 

The noise in a system can come from fluctuations in laser intensity, see ([5]), stray magnetic fields, fluctuations in 
currents producing the field and especially the ability to make a desired magnetic field in the lab. 
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A common approach for allowing noise in the system is to stretch the field ([7]) by taking a = rj • a max , where rj is 
a parameter (usually 0.5 < rj < 0.7), termed the 'Design Parameter'. This increases the length of the required field, 



-max — 2ria max ' 

Inserting in (jHJ yields, 

a(z) = (l-l) 7 , fl 7 ■ (9) 

For rj = 1, we have a = 0, not allowing for any noise in the system. For rj < 1; at the slower entrance, the allowed 
noise is (1 — rj) a ma x/vi, increasing monotonously with position along the slower. In most systems, and up to a good 
approximation, the noise is independent of position. Thus, from ©, the constant deceleration approach is not ideal. 

We now introduce two alternative approaches and show that they outperform again the traditional constant decel- 
eration approach. 

IV. ANALYTICAL APPROACH 

In this section we present a new approach for designing an optimal field profile for a system with position independent 
noise. 

For a (z) = a, the general solution to equation © reads, 



g , )= hka^ax , w 

fx' a 



*°—l]e 



ft,, 



1 , (10) 



where W is the Lambert W function defined as the function which solves z = W (z) e w ^ 2 
For a — > it is easy to check using a series expansion near — e -1 , 



W(z)a -1 + V2 (exz+1), (11 

that, 

hk 



B(z)^'jjyjv?-2za ma x, (12) 
which fits Q. 

As with the derivation of the traditional approach, we insert (0| in (TIT))) to obtain the relationship between length 
and maximal capture velocity, 

/ -, VfCX \ 

h-Vf 




(13) 

Q*max /J / 

In order to elucidate the advantages of such an approach let us consider a typical example. We consider the case of 
a metastable neon beam, interacting with a laser of intensity s = 3. The optical transition of interest, 3 P 2 ~~ ► 3 D 3 , 
has k = 9.8 • 10 6 to _1 , and V = 8.2 (2tt) Mhz. From |5]), a max = 6 • 10 5 m/s 2 . The magnetic moment of the transition 
is // = /is, where is the Bohr magneton. 

For a desired slower length length of z max = lm, a final velocity of, w/ = 0, and in the presence of noise that 
allows a design parameter of rj — 0.5, the maximal initial velocity one can slow using the conventional approach is 
Vi = ^2z rja max = 775 m/s. 

The minimum noise parameter is, a = (1 — rj) a ™ ax = 387s -1 . Substituting into (|14[) we obtain the maximum 
velocity allowed with the analytical approach, vi = 854 m/s, significantly larger. 

On the other hand, for a desired initial velocity of vi = 775 m/s, and rj — 0.5, the required length obtained from 
(IT5]) is z m ax — 0.77 m. Allowing a reduction of ^25% is the slower length. 

It is clear that the analytical approach yields better results, both in maximal capture velocity and in length. A 
Zeeman slower designed to have this field profile can be tailored for a higher capture velocity, less decollimation, or a 
combination of both. Fig. [T] shows the three different field profiles described in the text. 
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FIG. 1: (color online) A comparison of the three calculated Zeeman slower fields. The constant acceleration (square root) 
field (blue, dotted), the Lambert field designed for a known slower length (red, dashed), and the Lambert field designed for a 
specified capture velocity (green, solid). 



V. SYSTEM OPTIMIZATION APPROACH 



In the previous section we described an approach to slower design which allows an optimal field profile for a position 
independent noise profile. This fits most of the 'tapered solenoid' systems, where it is quite easy to make a desired 
target field and the noise is mostly from fluctuations in the current. 

In recent years, several systems have been designed, including our planned apparatus which have more degrees 
of freedom on the one hand, but are harder to fit to a target field on the other. We now show how to optimize a 
known system, using our own as an example. 

As a first step we write the magnetic field as a function of the controllable system parameters (currents, winding 
numbers, step-motor position (see [Hj) etc.), 

B = b(o,z\ (15) 

Now some magnetic field measurements should be done in order to determine how well this calculated (or simulated) 
field represents the real field in the lab. Any inconsistency between the two should be considered as a part of the 
'noise'. 

From (|8|), assuming a is position dependent, we can quantify a 'noise parameter' in units of G/m as a function of 
the controllable parameters, 

Y ^)- h i^tA B ^)) (16) 

We use the minimum noise parameter to avoid the atoms falling out of the adiabatic following condition, Eq. ([6]) , 
at any point along the slowing process. 

As a next step we measure or simulate fluctuations in the controllable parameters, in order to determine how they 
affect the 'noise parameter, and select the smallest, Yq, that takes them into account. This will be our optimization 
constraint, 

Y (e) > Y . (17) 
The yield function will be the parameter we wish to maximize. 
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Let us consider a specific system. Our group is building a midfield-zero modular Zeeman slower, consisting of several 
independently controlled, identical coils. We select a detuning of, So — 400 (2tt) Mhz, which from ^ determines the 
velocity in the zero crossing point to be, v c = So/k = 256to/s. Following (fl~5| we write, 



(18) 



where /„ is the current in the n th coil, and d = 77 mm is the distance between coil centers. The function / is an 
analytical fit to the field of a single coil at a current of 1 A. Comparing this calculated field to the one we measure at 
the lab gives good agreement up to ~ 0.5%. To be on the safe side we allow about 1% noise in the system. Varying 
the currents by up to 1% corresponds in our system to approximately Yq — 1 G /cm 

Once we have determined the constraint we can maximize the capture velocity. Our yield function is, 



— max 

hk 0<z<z mt 



B Lz 



So 
k ' 



(19) 



We are thus fitting for the maximum of a constrained (eq (|17p) nonlinear multivariable function (fl"§|) . We use Matlab's 
'fmincon' built in function which uses the 'active set' algorithm in a similar manner to the one described by [l3j . 
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FIG. 2: (color online) A comparison of the three fitted Zeeman slower fields for a constant length slower. The constant 
acceleration (square root) field (blue, dotted), the Lambert field designed for a known slower length (black, dashed), and the 
optimal field (red, solid). The coil center positions are also shown (orange, dots) 



Figure [5] shows the result of optimizing for a fixed slower length of 12 coils compared to the results of fitting to 
a square root target field and to a Lambert target field, all with with the same maximal noise Yq = 1 G j cm. The 
oscillations in the field profile stem from the discrete nature of the system and are taken into account implicitly 
in the noise parameter. In all cases the adiabatic following condition is satisfied. The optimized maximum field is 
506G, corresponding to a capture velocity of Vi — 710ra/s, whereas the Lambert field allows Mj = 682m/s, and the 
traditional (square root) field only permits Vi = 654 m/s. 

We may also use this method to design a shorter slower. We take for Vi — 654 m/s, which is the maximal velocity 
for a 12-coil slower designed with the traditional field approach to withstand noise of Yb = 1 G/cm. At this initial 
velocity we can use the Lambert target field with 9 coils, while still maintaining Yq>1 G/cm. When we optimize the 
system under the constraint, 



Vi < — max 

Hk 0<z<z ml 



B Lz 



So 
k ' 



(20) 
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and the yield function Y \ lj , we arrive at a slower with 8 coils. The results, along with the coil center positions, are 
shown in Fig. [3] 




z[mm] 

FIG. 3: (color online) A comparison of the three fitted Zeeman slower fields for a fixed capture velocity. The constant 
acceleration (square root) field (blue, dotted), the Lambert field (black, dashed), and the optimal field (red, solid). The coil 
center positions are also shown for each of the three fields (circles, triangles, and squares, respectively) 



VI. SUMMARY 

In summary we have introduced two new approaches for the design of a Zeeman slower which rely on optimal 
compliance with the adiabatic following condition. These approaches allow to optimize system parameters, such as 
the capture velocity and the slower length. 

We have shown that these two approaches yield better results than the traditionally used square-root (constant 
deceleration) approach. We note that the two approaches are complementary in that the first is useful for the design 
a system and the second one is useful for optimizing it after some measurements, such as the magnetic field and the 
noise, have been taken. 
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